Vector boson mass generation without new fields 
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Previously a model of only vector fields with a local U(l) (8)SU(2) symmetry was introduced for 
which one finds a massless U(l) photon and a massive SU(2) vector boson in the lattice regular- 
ization. Here it is shown that quantization of its classical continuum action leads to perturbative 
renormalization difficulties. But, non-perturbative Monte Carlo calculations favor the existence of 
a quantum continuum limit. 
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Electromagnetic currents plus charged and neutral 
weak currents with the corresponding vector bosons 
{■-f,W^,Z^), are needed for a theory that accommo- 
dates simultaneously weak parity violation and electro- 
magnetic parity conservation Explicit breaking of 
gauge symmetry through massive vector bosons can be 
avoided by the Higgs mechanism which leads to de- 
sired features such as perturbative renormalizability Q. 

Nevertheless, the introduction of the Higgs particle 
into a theory in which all matter fields are fermions with 
interactions mediated by vector bosons remains quite ad- 
hoc and the quadratic divergence of its self-energy causes 
a fine tuning problem This provides an opening for 
proposing new physics solutions of which supersymmetry 
and string theory are most popular. Though new physics 
should emerge on the way to the Planck scale, there ap- 
pear no strong reasons to expect experimental signals for 
it on the LHC energy scale. Occam's razor suggests to 
stay with fermions and vector bosons. 

Notably, the original arguments from the 1960s and 
early 1970s rely on perturbation theory. Our non- 
perturbative understanding of quantum field theory 
(QFT) developed later. A milestone was Wilson's [5] 
formulation of lattice gauge theory (LGT) in 1974 and 
LGT Monte Carlo (MC) calculations started in earnest 
during the early 1980s after pioneering papers by Creutz 
and others 0] . In recent years that the author suggested 
to address the problem of vector boson mass generation 
from scratch within the non-perturbative LGT frame- 
work 7]. This appears to be of theoretical interest, in- 
dependently of the question whether such a scenario is 
eventually realized by nature or not. 

Using Wilson's regularization, = ex.'p{igbaB ^) , 
— f • 6^/2, where r^, i — 1,2,3 are the Pauli ma- 
trices, the SU(2) lattice action reads 

p 

where the sum is over all plaquettes and Vp are oriented 
products of SU(2) matrices around the plaquette loop. 
E.g., for a plaquette in the fiv, 1^ 7^ plane Vp becomes 



(a lattice spacing, x — na with n an integer 4- vector): 

V(a^) - i:r[V^,{x)V,{x + fia)Vl{x + i'a)V}{x)\ . 

Due to non-perturbative effects there are no massless par- 
ticles in the spectrum of SU(2) LGT. The self interaction 
of the gauge fields creates a spectrum of massive glueballs 
and one of them may be used to set the mass scale. Cou- 
pled fermions are confined, while the leptons are found 
as free particles. So, we need SU(2) LGT in a deconfined 
phase. This can be achieved by increasing the physical 
temperature 0] , but in a fundamental theory we have to 
stay at zero temperature. 

The order parameter for the deconfining phase tran- 
sition is the expectation value of the Polyakov loop. 
Polyakov loops are traces of products of SU(2) matrices 
along straight lines closed by periodic boundary condi- 
tions. On a finite lattice one finds from confined to de- 
confined a transition between a single peak and a double 
peak distribution. A coupling of parallel SU(2) matrices, 

ReTr [V,{x + (ia)V^{x)\ , 

is well suited to align Polyakov loops, but breaks SU(2) 
gauge invariance. To rescue local U(l) ® SU(2) invari- 
ance of matter fields a concept of extended gauge in- 
variance was introduced 0, which requires to introduce 
additional vector fields. 

Defining the U(l) field as 2 x 2 matrix Up, — 
expiigaaAfj), A^^ = tq ap/2 (tq unit matrix), we consider 
in the forthcoming the vector field lattice action 

S^^J2^eT,Up + ^J2T^rVp + ^J2sf.'. W 

p p 

where the third sum includes identical j.i = u indices and 

Sfj^ =ReTr [Up{x)V,,{x + fm)Ul{x + va)V}{x)\ . 

The relations Pa — l/(?a and /S^ — '^j g\ define and 
/9b of (U) through the bare couplings constants and A is 
a new free parameter. Properties as function of A have 
been investigated by MC calculations J\ . After fixing /3a 
in the Coulomb phase of U(l) LGT and /3b in the scaling 
region of confined SU(2) LGT, one finds for small A the 
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same results as for A = 0: A massless U(l) photon and 
confined SU(2) gauge theory with a glueball spectrum. 
Increasing A, a strong first order phase transition takes 
place. The U(l) photon survives the transition mass- 
less, while SU(2) is then in a deconfined phase with a 
massive vector boson triplet. Central questions are then 
about 1. Perturbative Renormalizability and 2. Existence 
of a Quantum Continuum Limit. Both remained beyond 
the scope of 0. New insights are reported here after 
discussing a novel way to ensure local U(l) ® SU(2) in- 
variance. 

Let us consider U(l) (g) SU(2) field configurations 
{U'^{na)} for which the SU(2) part is a gauge transforma- 
tion of the zero field B^j_{na) = and {V^{na)} for which 
the U(l) part is a gauge transformation of the zero field 
Ai^{na) = 0. With G G U(l) ® SU(2) each set is mapped 
onto itself by 

C//,(na) ^ G{na)Ul{na)G-^{na + jia), (2) 
V^(na) ^ G{na)Vl^{na)G-^{na + fia), (3) 

which we call extended gauge transformations. Replacing 
and by their primed versions, the action ([T]) is in- 
variant under these transformations. S^'^ is the simplest 
example of Wilson loops that mix U'^ and V"^ matrices, 
which are now invariant operators. 
The partition function of the model is 

^ = (4) 

n /I— 1 

where the integrations are over {U'^{na)} and {V^{na)} 
defined above. They are easily implemented in a MC 
calculation. Let us use the notation proper parts for 
the U(l) factor of the U'^ and the SU(2) factor of the 
matrices and the notation gauge parts for the SU(2) 
factor of the U'^ and the U(l) factor of the V"^ matri- 
ces. The proper parts can be updated in the usual way. 
Specifically, a biased Metropolis-heatbath algorithm [9] 
was used in the simulations. Updates of the gauge parts 
transform all matrices emerging at a site n according to 

C/;(na) ^ G^ina) U'^{na) , ^/(n) ^ Gi(na) V'^{na) , (5) 

and all matrices on links ending at n according to 

Ufj,{na^fla) — >■ U ^{na — jla) G2^ ina) ^ 

V^{na — (la) Vi^{na — fia) G^^ (na) , (6) 

where G2 and Gi are, respectively, SU(2) and U(l) matri- 
ces drawn with the group measure. These updates change 
[lo| the action ([T]), so that a Metropolis algorithm will 
have an acceptance rate in the range (0, 1]. 

Changes of the action under ([5]) can be undone by ap- 
propriate updates of the proper parts of the matrices. 
Therefore, we can calculate operators that are invari- 
ant under extended gauge transformations in any fixed 



gauge, i.e., omitting updates of the form ([5]). For MC 
calculations it is convenient to assign zero fields to the 
gauge parts, which we call proper gauge. The purpose 
of extended gauge invariance is to allow an initial La- 
grangian with local U(l) ® SU(2) invariance of matter 
fields without explicit breaking by a vector boson mass 
term. Including matter fields our Lagrangian in the clas- 
sical continuum limit is (using Euclidean notation) 

L^Jj {i^^Dl^ - m) {i^^Dl - m) ^ , (7) 

--Tr ( F'^ F"^ ) — - ( F'' F'' ) - — Tr ( F^^'^ padd\ 

Here D"^ ^ df_, + igaA'^ and D^^^ + igbB'^ are gauge 
covariant derivatives and the additional field tensor is 

- 9bdf.Bl - gad.A'^ + I gag, [A;,, B',\ . (8) 

The fermion field '0 is assumed to be a doublet and 
the Lagrangian is invariant under the local U(l) ® SU(2) 
symmetry transformations tJj ^ G ip with the contin- 
uum limit of extended gauge transformations for the 
vector fields being A'^ GA'^G'^ + i{d^G)G-^ / ga, 
B'^ GB'^G-^ + i{df,G)G-^/gb 7]. This is the rea- 
son for the occurrence of two V' terms in the Lagrangian. 

1. In the proper gauge, S^, — > A^, one gets 

L^'' = id,a.f - ^ {d,Kf ■ (9) 

These pieces are found in U(l) and SU(2) effective La- 
grangians, which are usually obtained by integrating 
over gauge transformations with a Gaussian weighting 
function. With the identifications ^ = 8/(A5^) and 
^ = 8/(Ag2), respectively, Eqn. (9.56) and (16.34) of [ul. 
However, (16.43) comes here without the Faddeev-Popov 
ghost fields. Therefore [ll[, the tree approximation 
is non-unitarity because of transitions to longitudinal 
modes, which require massive vector bosons while there 
is no explicit mass term in the Lagrangian (O, while the 
lattice regularization is unitary to the extent that one 
can prove reflection positivity. On the 1-loop level the 
vector boson self-energy is divergent, generating an infi- 
nite mass. These properties render the model ill-defined 
in conventional perturbation theory. 

2. One jnay expect that the vector boson mass amw 
found in [7] is also non-perturbatively divergent. Then, 
the lattice regularization would not allow for a quantum 
continuum limit amw 0. Instead, amw has to stay fi- 
nite amw > ammin > in a smoothly connected range of 
couplings, eventually bounded by first order phase tran- 
sitions. We investigate here the line 

13a = \, /3b = 2A, A ^00 (10) 

for which one could envision an approach to a quantum 
continuum limit in analogy to the behavior of asymp- 
totically free non-Abelian gauge theories. The result is 
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TABLE I: Mass estimates on lattices and infinite vol- 

ume extrapolations according to Eq. (|13p . 



A 


A^min 


amw{^, 14) 


amiv(A, oo) 


Q 


1.1 


6 


0.2659 (10) 


0.2658 (10) 


0.13 


4.0 


4 


0.1245 (10) 


0.1180 (16) 


0.31 


8.0 


4 


0.0905 (12) 


0.0876 (15) 


0.44 


12.0 


4 


0.0740 (11) 


0.0719 (14) 


0.31 


16.0 


4 


0.0675 (14) 


0.0653 (13) 


0.36 


20.0 


4 


0.0597 (14) 


0.0552 (16) 


0.42 


24.0 


4 


0.0523 (10) 


0.0500 (15) 


0.11 


28.0 


6 


0.0519 (10) 


0.0487 (18) 


0.26 


32.0 


4 


0.0480 (13) 


0.0448 (18) 


0.96 



that our fits to the scenarios amw arrinun > versus 
amw — >■ prefer the latter. 

Our mass spectrum calculations were performed on lat- 
tices of size N^Nt, Nt ^ N. For each value of A we have 
first to extrapolate the infinite volume limit N ^ oo 
amwi^,N), denoted by amwW — amwi^, oo). Subse- 
quently, we fit amwW so that a A — oo extrapolation 
along the line defined by ((TO)) can be performed. 

Our masses are deduced from correlation functions c{t) 
of suitable trial operators by performing the usual two 
parameter cosh fits 

c{t) = ai [exp{—amwt)+e:cp{~amw{Nt~t))] (11) 

for a range of integers < ti < t < t2- Here the trial 
operators 

VK,,^(x) = -^Tr [nW^^ix)] , (12) 

are employed, where in slight deviation from 7] a U(l) 
phase is included, Wf^(x) = Uj^{x) Vf^{x) (no summation 
over /i). As the previously used operator, it becomes 
gauge invariant in combination with (static) fermion or 
boson fields, compare (6.20) of il2i] . 

In addition correlations between trial operators for the 
U(l) photon and SU(2) glueball masses in the plaquette 
representations of the cubic group were calculated. Es- 
timates of the U(l) photon mass are for all A consistent 
with zero, while there are no convincing signals in the 
glueball channels. This is similar to the results reported 
in j?! for one choice of coupling constant values. 

Simulations were carried out on lattices of size N — 
4, 6, 8, 10, 12, 14 with Nt in the range [48,96] and A val- 
ues as given in table HI For each A the 3-paramcter fit 

amw{^,N) — amw{^,oo)+a2{X) exp{~a3{X)N) (13) 

was performed to derive an infinite volume estimate 
amw{X, oo). These fits are shown in Fig. [TJ The ex- 
trapolations are collected in table |T] together with the 
goodness Q of each fit and the estimates on our largest 



0.4 ^ 




4 6 8 10 12 14 

N 



FIG. 1: Fits of amivi^, N). The A values correspond in up 
to down order to the curves. 
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FIG. 2: Correlation functions for the A = 20 mass estimates. 
The A'^ values correspond in up to down order to the curves. 



14'^iV( lattices. Error bars are given in parenthesis and 
refer to the last digits in the number before. In two cases 
data from the smallest 4'^A''f lattice were omitted from 
the fit for consistency reasons. We indicate with Nn^n 
the size of the smallest lattice included in the fit. 

To give an example of the numerical quality of the 
correlation functions (ITT]) we depict them in Fig. [2] for 
our lattices at A = 20. In stark contrast to the noise 
one encounters for glueball correlations in pure lattice 
gauge theories, these are beautiful strong correlations. 
One can easily follow them over more than 30 lattice 
spacing, though the estimates for our 14^A't lattices are 
already rather time consuming. Relying on a statistics 
of 500,000 sweeps, they run with the present single pro- 
cessor code one week on an Intel 17 CPU. 

We are now prepared to discuss the A — > oo behavior 
of the amvi/(A) = amwiX, oo) values of table HI A num- 
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TABLE II: Fits of the amw(A) = amw{^, oo) values of ta- 
ble m The first column gives the number of parameters, the 
last column the goodness of the fit. 



par # 


function 


amvK(oo) 


Q 


2 


/i(A) = ai exp(-a2A) 








3 


/2(A) = ao + /i(A) 


0.05783 (66) 





3 


-f / \ \ \ — CIO ( \\ 

/3(A) = aiA 2 exp(-a3A) 





1.2 X 10 


4 


/4(A) = ao + /3(A) 


0.03 (18) 


1.1 X 10"" 


4 


/5(A) = /3(A) (l + a4/A) 





0.26 


5 


/6(A) = ao + /4(A) 


0.02 (11) 


0.20 


3 


/7(A) = aiA-''^(l + a4/A) 





0.036 


4 


/8(A) = ao + /7(A) 


-0.1 (1.1) 


0.065 


2 


/9(A) = aiA-"^ 





4.7 X 10"^^ 


10 


/10(A) = ao + /9(A) 


0.0241 (26) 


2.0 X 10"'' 



s 

E 

CO 



0.3 



Mass estimates 
Fit 



Goodness of fit Q=0.26 




value of the factor 03 in the exponent is small, one 
wonders whether a power law alone suffices to describe 
mn/(A). The functional forms /7(A) to /10(A) test this. 
While the last two of them are bad, this is less obvious 
for /7(A) and fsW- Recall, under the assumption that 
the form of a fit is correct, Q is the probability for the 
discrepancy between the fit and the data. 

The presented MC calculations indicate divergence of 
the correlation length ^/a — ^ 00, ^ = for A — 

00, contradicting perturbation theory and supporting a 
quantum continuum limit. Of course, we cannot exclude 
that for larger systems and coupling constants the behav- 



ior may turn around and support mvi/(A) 



> 0. 



While this is correct, the perturbative approach is in 
essence vulnerable to similar criticism. Though it served 
us well, there is no proof that perturbation theory de- 
scribes the true nature of a QFT. 

Besides moving on to new horizons, we should perhaps 
keep an open mind for the third logical possibility that a 
deeper understanding of conventional QFT could unveil 
new models with massive vector bosons. The unexpected 
numerical results of this paper provide no answers, but 
indicate that it is worthwhile to continue this line of work. 
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Some of the computer programs used rely on collabo- 
rations with Alexei Bazavov. 



FIG. 3: Best fit of amw(A). 



ber of fits, either enforcing aTOn/(A) for A — > cxd or 
allowing for a free parameter amw {00) were tried and 
are compiled in table HH The fit forms /i(A) to /4(A) are 
in disagreement with the data as signaled by very small 
Q values in the last column. A zero means that Q is so 
small that the precise value cannot be tracked within our 
rounding errors of about 10"^°. Most convincing is our 
4-parameter fit 

amwW = /5(A) = aiA"^ exp(a3A) (1 a4/A) (14) 

with a well acceptable goodness of fit Q = 0.26. The 
values of its parameters are oi = 0.111(47), 02 — 

-0.16(18), 03 = -0.0139(75), 04 = 1.6 (^^g e)" ^^^^ 

adding amw {00) as 5th parameter ao, the function /6(A) 
is obtained for which the goodness of fit goes slightly 
down, indicating that we are overfitting, and the esti- 
mated ao = amw (00) includes zero. 

The functional form has similarities with the 

asymptotic freedom behavior in LGT. As the estimated 
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